Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score--the CoLaus Study.

AIMS/HYPOTHESIS: Several susceptibility genes for type 2 diabetes have been discovered recently. Individually, these genes increase the disease risk only minimally. The goals of the present study were to determine, at the population level, the risk of diabetes in individuals who carry risk alleles within several susceptibility genes for the disease and the added value of this genetic information over the clinical predictors. METHODS: We constructed an additive genetic score using the most replicated single-nucleotide polymorphisms (SNPs) within 15 type 2 diabetes-susceptibility genes, weighting each SNP with its reported effect. We tested this score in the extensively phenotyped population-based cross-sectional CoLaus Study in Lausanne, Switzerland (n = 5,360), involving 356 diabetic individuals. RESULTS: The clinical predictors of prevalent diabetes were age, BMI, family history of diabetes, WHR, and triacylglycerol/HDL-cholesterol ratio. After adjustment for these variables, the risk of diabetes was 2.7 (95% CI 1.8-4.0, p = 0.000006) for individuals with a genetic score within the top quintile, compared with the bottom quintile. Adding the genetic score to the clinical covariates improved the area under the receiver operating characteristic curve slightly (from 0.86 to 0.87), yet significantly (p = 0.002). BMI was similar in these two extreme quintiles. CONCLUSIONS/INTERPRETATION: In this population, a simple weighted 15 SNP-based genetic score provides additional information over clinical predictors of prevalent diabetes. At this stage, however, the clinical benefit of this genetic information is limited.

Effect of single nucleotide polymorphisms in cytochrome P450 isoenzyme and N-acetyltransferase 2 genes on the metabolism of artemisinin-based combination therapies in malaria patients from Cambodia and Tanzania.

The pharmacogenetics of antimalarial agents are poorly known, although the application of pharmacogenetics might be critical in optimizing treatment. This population pharmacokinetic-pharmacogenetic study aimed at assessing the effects of single nucleotide polymorphisms (SNPs) in cytochrome P450 isoenzyme genes (CYP, namely, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, and CYP3A5) and the N-acetyltransferase 2 gene (NAT2) on the pharmacokinetics of artemisinin-based combination therapies in 150 Tanzanian patients treated with artemether-lumefantrine, 64 Cambodian patients treated with artesunate-mefloquine, and 61 Cambodian patients treated with dihydroartemisinin-piperaquine. The frequency of SNPs varied with the enzyme and the population. Higher frequencies of mutant alleles were found in Cambodians than Tanzanians for CYP2C9*3, CYP2D6*10 (100C → T), CYP3A5*3, NAT2*6, and NAT2*7. In contrast, higher frequencies of mutant alleles were found in Tanzanians for CYP2D6*17 (1023C → T and 2850C → T), CYP3A4*1B, NAT2*5, and NAT2*14. For 8 SNPs, no significant differences in frequencies were observed. In the genetic-based population pharmacokinetic analyses, none of the SNPs improved model fit. This suggests that pharmacogenetic data need not be included in appropriate first-line treatments with the current artemisinin derivatives and quinolines for uncomplicated malaria in specific populations. However, it cannot be ruled out that our results represent isolated findings, and therefore more studies in different populations, ideally with the same artemisinin-based combination therapies, are needed to evaluate the influence of pharmacogenetic factors on the clearance of antimalarials.

Coevolution of restriction factors and retroviruses

Les virus exploitent la machinerie cellulaire de l'hôte pour se répliquer. Ils doivent s'adapter pour infecter la cellule hôte de manière optimale tout en échappant à la vigilance du système de défense de l'hôte. Ainsi l'hôte et les virus se livrent à de constantes batailles évolutives. Mon travail de thèse a porté sur l'étude des signatures évolutives de facteurs de l'hôte agissant comme des 'facteurs de restriction' en bloquant la réplication rétrovirale chez les primates. Plus spécifiquement, mon travail a visé à utiliser des données évolutives pour renseigner les analyses fonctionnelles et la biologie. Nous avons étudié le facteur anti-VIH-1 nommé TRIM5a (i) chez les prosimiens pour mieux comprendre son rôle dans le contrôle d'un lentivirus endogène, (ii) dans son activité contre d'autres anciennes infections représentées par des rétrovirus endogènes humains et (iii) en tant que protéine capable de générer des mutants de la capside. Premièrement nous nous sommes intéressés à TRIM5a chez deux espèces de lémuriens dont Microcebus murinus qui porte le lentivirus endogène PSIV dans son génome depuis plusieurs millions d'années,. Nous avons observé que TRIM5a chez M. murinus a un spectre d'activité antivirale réduit à l'opposé de TRIM5a chez le Lemur catta - non porteur du PSIV endogène - qui bloque une large variété de rétrovirus dont le PSIV. De ce fait TRIM5a aurait pu contribuer à protéger certaines espèces de lémuriens vis-à-vis d'anciennes infections par le PSIV. A l'inverse du PSIV, des virus dérivés des rétrovirus endogènes humains HERV-K and HERV-H se sont révélés largement résistants à l'inhibition par TRIM5a. Ces données illustrent une absence de protection par TRIM5a face à d'autres anciennes infections rétrovirales. Puis, pour évaluer l'impact de la protéine TRIM5a humaine sur le VIH-1, nous avons testé l'effet de mutations des résidues sous sélection positive dans la capside du VIH-1 sur l'inhibition par TRIM5a. Nos résultats montrent que TRIM5a ne jouerait pas un rôle significatif dans l'évolution de la capside du VIH-1. Enfin notre travail a porté sur le facteur anti-VIH-1 SAMHD1 récemment découvert, que nous avons séquencé chez 25 espèces de primates. L'analyse évolutive des sites sous sélection positive et des expériences fonctionnelles ont permis d'identifier le domaine de SAMHD1 interagissant avec la protéine lentivirale Vpx. De même que d'autres protéines virales contrecarrent les facteurs de restriction en les menant à la dégradation, nous avons observé que Vpx induit la dégradation de SAMHD1 de manière spécifique à l'espèce. Ces découvertes contribuent à comprendre comment les facteurs de restriction et les virus co-évoluent pour se neutraliser l'un l'autre. - Viruses hijack the host cellular machinery to replicate. They adapt to infect optimally host cells while escaping host defense systems. Viruses and the host coevolve in an evolutionary struggle. My thesis work has been devoted to study the evolutionary signatures of host factors acting as restriction factors that block retroviral replication in primates. Specifically, my work aimed at using evolutionary data to inform functional analyses and biology. We studied the anti-HIV-1 factor TRIM5a (i) in prosimians to better understand its possible role in the control of an endogenous lentivirus, (ii) in its activity against other ancient infections - as represented by HERVs, and (iii) as a protein capable of generating escape mutants in the viral capsid. First, my work focused on two lemur species, one of which was the gray mouse lemur that carries the endogenous lentivirus PSIV integrated in its genome for several million years. TRIM5a from gray mouse lemur exhibited a limited antiviral spectrum as opposed to TRIM5a from ring-tailed lemur - not a host of PSIV - that is able to block diverse retroviruses notably PSIV. These results support the possible contribution of TRIM5a in protecting lemur species from ancient infection by PSIV. In contrast, chimeric viruses derived from two human endogenous retroviruses were broadly resistant to TRIM5a-mediated restriction, suggesting TRIM5a lack of activity against other types of ancient infections. To evaluate the recent impact of human TRIM5a on HIV-1 evolution, we tested whether variants at positively selected sites in the HIV-1 capsid affected the ability of human TRIM5a alleles to restrict HIV-1. Our results indicate that TRIM5a does not play a significant role in the evolution of HIV1 capsid. At last, our work concentrated on the newly discovered anti-HIV-1 restriction factor SAMHD1. We determined its coding sequence in a panel of 25 species of primates. Evolutionary analyses of positively selected sites in SAMHD1 domains and functional assays identified the domain of SAMHD1 interacting with the lentiviral protein Vpx. Similar to other viral countermeasures targeting cellular restriction factors, Vpx was responsible of the degradation of SAMHD1 orthologs in a species-specific manner. These findings contributed to understanding how restriction factors and viruses evolve to counteract each other.

The protein phosphatase 7 regulates phytochrome signaling in Arabidopsis.

The psi2 mutant of Arabidopsis displays amplification of the responses controlled by the red/far red light photoreceptors phytochrome A (phyA) and phytochrome B (phyB) but no apparent defect in blue light perception. We found that loss-of-function alleles of the protein phosphatase 7 (AtPP7) are responsible for the light hypersensitivity in psi2 demonstrating that AtPP7 controls the levels of phytochrome signaling. Plants expressing reduced levels of AtPP7 mRNA display reduced blue-light induced cryptochrome signaling but no noticeable deficiency in phytochrome signaling. Our genetic analysis suggests that phytochrome signaling is enhanced in the AtPP7 loss of function alleles, including in blue light, which masks the reduced cryptochrome signaling. AtPP7 has been found to interact both in yeast and in planta assays with nucleotide-diphosphate kinase 2 (NDPK2), a positive regulator of phytochrome signals. Analysis of ndpk2-psi2 double mutants suggests that NDPK2 plays a critical role in the AtPP7 regulation of the phytochrome pathway and identifies NDPK2 as an upstream element involved in the modulation of the salicylic acid (SA)-dependent defense pathway by light. Thus, cryptochrome- and phytochrome-specific light signals synchronously control their relative contribution to the regulation of plant development. Interestingly, PP7 and NDPK are also components of animal light signaling systems.

Molecular variation at a candidate gene implicated in the regulation of fire ant social behavior

The fire ant Solenopsis invicta and its close relatives display an important social polymorphism involving differences in colony queen number. Colonies are headed by either a single reproductive queen (monogyne form) or multiple queens (polygyne form). This variation in social organization is associated with variation at the gene Gp-9, with monogyne colonies harboring only B-like allelic variants and polygyne colonies always containing b-like variants as well. We describe naturally occurring variation at Gp-9 in fire ants based on 185 full-length sequences, 136 of which were obtained from S. invicta collected over much of its native range. While there is little overall differentiation between most of the numerous alleles observed, a surprising amount is found in the coding regions of the gene, with such substitutions usually causing amino acid replacements. This elevated coding-region variation may result from a lack of negative selection acting to constrain amino acid replacements over much of the protein, different mutation rates or biases in coding and non-coding sequences, negative selection acting with greater strength on non-coding than coding regions, and/or positive selection acting on the protein. Formal selection analyses provide evidence that the latter force played an important role in the basal b-like lineages coincident with the emergence of polygyny. While our data set reveals considerable paraphyly and polyphyly of S. invicta sequences with respect to those of other fire ant species, the b-like alleles of the socially polymorphic species are monophyletic. An expanded analysis of colonies containing alleles of this clade confirmed the invariant link between their presence and expression of polygyny. Finally, our discovery of several unique alleles bearing various combinations of b-like and B-like codons allows us to conclude that no single b-like residue is completely predictive of polygyne behavior and, thus, potentially causally involved in its expression. Rather, all three typical b-like residues appear to be necessary.

Direct binding of peptides to MHC class I molecules on living cells. Analysis at the single cell level.

To directly assess the binding of exogenous peptides to cell surface-associated MHC class I molecules at the single cell level, we examined the possibility of combining the use of biotinylated peptide derivatives with an immunofluorescence detection system based on flow cytometry. Various biotinylated derivatives of the adenovirus 5 early region 1A peptide 234-243, an antigenic peptide recognized by CTL in the context of H-2Db, were first screened in functional assays for their ability to bind efficiently to Db molecules on living cells. Suitable peptide derivatives were then tested for their ability to generate positive fluorescence signals upon addition of phycoerythrin-labeled streptavidin to peptide derivative-bearing cells. Strong fluorescent staining of Db-expressing cells was achieved after incubation with a peptide derivative containing a biotin group at the C-terminus. Competition experiments using the unmodified parental peptide as well as unrelated peptides known to bind to Kd, Kb, or Db, respectively, established that binding of the biotinylated peptide to living cells was Db-specific. By using Con A blasts derived from different H-2 congenic mouse strains, it could be shown that the biotinylated peptide bound only to Db among > 20 class I alleles tested. Moreover, binding of the biotinylated peptide to cells expressing the Dbm13 and Dbm14 mutant molecules was drastically reduced compared to Db. Binding of the biotinylated peptide to freshly isolated Db+ cells was readily detectable, allowing direct assessment of the relative amount of peptide bound to distinct lymphocyte subpopulations by three-color flow cytometry. While minor differences between peripheral T and B cells could be documented, thymocytes were found to differ widely in their peptide binding activity. In all cases, these differences correlated positively with the differential expression of Db at the cell surface. Finally, kinetic studies at different temperatures strongly suggested that the biotinylated peptide first associated with Db molecules available constitutively at the cell surface and then with newly arrived Db molecules.

Immunochip SNP array identifies novel genetic variants conferring susceptibility to candidaemia.

Candidaemia is the fourth most common cause of bloodstream infection, with a high mortality rate of up to 40%. Identification of host genetic factors that confer susceptibility to candidaemia may aid in designing adjunctive immunotherapeutic strategies. Here we hypothesize that variation in immune genes may predispose to candidaemia. We analyse 118,989 single-nucleotide polymorphisms (SNPs) across 186 loci known to be associated with immune-mediated diseases in the largest candidaemia cohort to date of 217 patients of European ancestry and a group of 11,920 controls. We validate the significant associations by comparison with a disease-matched control group. We observe significant association between candidaemia and SNPs in the CD58 (P = 1.97 × 10(-11); odds ratio (OR) = 4.68), LCE4A-C1orf68 (P = 1.98 × 10(-10); OR = 4.25) and TAGAP (P = 1.84 × 10(-8); OR = 2.96) loci. Individuals carrying two or more risk alleles have an increased risk for candidaemia of 19.4-fold compared with individuals carrying no risk allele. We identify three novel genetic risk factors for candidaemia, which we subsequently validate for their role in antifungal host defence.

DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe.

STUDY OBJECTIVE: Prior research has identified five common genetic variants associated with narcolepsy with cataplexy in Caucasian patients. To replicate and/or extend these findings, we have tested HLA-DQB1, the previously identified 5 variants, and 10 other potential variants in a large European sample of narcolepsy with cataplexy subjects. DESIGN: Retrospective case-control study. SETTING: A recent study showed that over 76% of significant genome-wide association variants lie within DNase I hypersensitive sites (DHSs). From our previous GWAS, we identified 30 single nucleotide polymorphisms (SNPs) with P < 10(-4) mapping to DHSs. Ten SNPs tagging these sites, HLADQB1, and all previously reported SNPs significantly associated with narcolepsy were tested for replication. PATIENTS AND PARTICIPANTS: For GWAS, 1,261 narcolepsy patients and 1,422 HLA-DQB1*06:02-matched controls were included. For HLA study, 1,218 patients and 3,541 controls were included. MEASUREMENTS AND RESULTS: None of the top variants within DHSs were replicated. Out of the five previously reported SNPs, only rs2858884 within the HLA region (P < 2x10(-9)) and rs1154155 within the TRA locus (P < 2x10(-8)) replicated. DQB1 typing confirmed that DQB1*06:02 confers an extraordinary risk (odds ratio 251). Four protective alleles (DQB1*06:03, odds ratio 0.17, DQB1*05:01, odds ratio 0.56, DQB1*06:09 odds ratio 0.21, DQB1*02 odds ratio 0.76) were also identified. CONCLUSION: An overwhelming portion of genetic risk for narcolepsy with cataplexy is found at DQB1 locus. Since DQB1*06:02 positive subjects are at 251-fold increase in risk for narcolepsy, and all recent cases of narcolepsy after H1N1 vaccination are positive for this allele, DQB1 genotyping may be relevant to public health policy.

Mutation screening of the glutamate cysteine ligase modifier (GCLM) gene in patients with schizophrenia.

BACKGROUND: Experimental evidences show that glutathione and its rate-limiting synthesizing enzyme, the glutamate-cysteine ligase (GCL), are involved in the pathogenesis of schizophrenia. Furthermore, genetic association has been previously reported between two single nucleotide polymorphisms lying in noncoding regions of glutamate cysteine ligase modifier (GCLM) gene, which specifies for the modifier subunit of GCL and schizophrenia. OBJECTIVE: We wanted to investigate the presence of GCLM true functional mutations, likely in linkage disequilibrium with the previously identified single nucleotide polymorphism alleles, in the same set of cases that allowed the detection of the original association signal. METHODS: We screened all the coding regions of GCLM and their intronic flanking vicinities in 353 patients with schizophrenia by direct DNA sequencing. RESULTS: Ten sequence variations were identified, five of which were not previously described. None of these DNA changes was within the GCLM coding sequence and in-silico analysis failed to indicate functional impairment induced by these variations. Furthermore, screening of normal controls and downstream statistical analyses revealed no significant relationship of any of these DNA variants with schizophrenia. CONCLUSION: It is unlikely that functional mutations in the GCLM gene could play a major role in genetic predisposition to schizophrenia and further studies will be required to assess its etiological function in the disease.

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10(-4)). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10(-3), n = 22,044), increased triglycerides (p = 2.6×10(-14), n = 93,440), increased waist-to-hip ratio (p = 1.8×10(-5), n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10(-3), n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10(-13), n = 96,748) and decreased BMI (p = 1.4×10(-4), n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.

Genome-wide association analysis identifies 20 loci that influence adult height.

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.

Genetic variation near IRS1 associates with reduced adiposity and an impaired metabolic profile.

Genome-wide association studies have identified 32 loci influencing body mass index, but this measure does not distinguish lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ∼2.5 million SNPs and body fat percentage from 36,626 individuals and followed up the 14 most significant (P < 10(-6)) independent loci in 39,576 individuals. We confirmed a previously established adiposity locus in FTO (P = 3 × 10(-26)) and identified two new loci associated with body fat percentage, one near IRS1 (P = 4 × 10(-11)) and one near SPRY2 (P = 3 × 10(-8)). Both loci contain genes with potential links to adipocyte physiology. Notably, the body-fat-decreasing allele near IRS1 is associated with decreased IRS1 expression and with an impaired metabolic profile, including an increased visceral to subcutaneous fat ratio, insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.

Strain-transcending Fc-dependent killing of Plasmodium falciparum by merozoite surface protein 2 allele-specific human antibodies.

It is widely accepted that antibody responses against the human parasitic pathogen Plasmodium falciparum protect the host from the rigors of severe malaria and death. However, there is a continuing need for the development of in vitro correlate assays of immune protection. To this end, the capacity of human monoclonal and polyclonal antibodies in eliciting phagocytosis and parasite growth inhibition via Fcγ receptor-dependent mechanisms was explored. In examining the extent to which sequence diversity in merozoite surface protein 2 (MSP2) results in the evasion of antibody responses, an unexpectedly high level of heterologous function was measured for allele-specific human antibodies. The dependence on Fcγ receptors for opsonic phagocytosis and monocyte-mediated antibody-dependent parasite inhibition was demonstrated by the mutation of the Fc domain of monoclonal antibodies against both MSP2 and a novel vaccine candidate, peptide 27 from the gene PFF0165c. The described flow cytometry-based functional assays are expected to be useful for assessing immunity in naturally infected and vaccinated individuals and for prioritizing among blood-stage antigens for inclusion in blood-stage vaccines.